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Abstract 

Following the discovery of two new scalar mesons, /o(1370) and /o(1500) at the 
Low Energy Antiproton Ring at CERN, we argue that the observed properties of 
this pair are incompatible with them both being QQ mesons. We show instead that 
/o(1500) is compatible with the ground state glueball expected around 1500 MeV 
mixed with the nearby states of the ++ QQ nonet. Tests of this hypothesis include 
the prediction of a further scalar state, /^(lSOO — 1800) which couples strongly to 
KK, r/r/ and nrf. Signatures for a possible tensor glueball at ~2 GeV are also 
considered. 
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1 Introduction 



Glueballs are a missing link of the standard model. Whereas the gluon degrees of freedom 
expressed in Lqcd have been established beyond doubt in high momentum data, their 
dynamics in the strongly interacting limit epitomised by hadron spectroscopy are quite 
obscure. This may be about to change as a family of candidates for gluonic hadrons 
(glueballs and hybrids) is now emerging [JJ, |2J. In this paper we shall argue that scalar 
mesons around 1.5 GeV, in particular the detailed phenomenology of / (1500) and its 
partner / (1370), suggest that a glueball exists in this region, probably mixed with nearby 
isoscalar members of the scalar nonet. This hypothesis may be tested in forthcoming 
experiments. 

In advance of the most recent data, theoretical arguments suggested that there may 
be gluonic activity manifested in the 1.5 GeV mass region. Lattice QCD is the best 
simulation of theory and predicts the lightest "primitive" (ie quenched approximation) 
glueball to be ++ with mass 1.55 ± 0.05 GeV ||. Recent lattice computations place 
the glueball slightly higher in mass at 1.74 ± 0.07 GeV || with an optimised value for 
phenomenology proposed by Teper|5[] of 1.57±0.09 GeV. That lattice QCD computations 
of the scalar glueball mass are now concerned with such fine details represents considerable 
advance in this field. Whatever the final concensus may be, these results suggest that 
scalar mesons in the 1.5 GeV region merit special attention. Complementing this has 
been the growing realisation that there are now too many ++ mesons confirmed for them 
all to be QQ states |, |, 0. 

At ~ 1.5 GeV there is a clear ++ signal, / (1500), in several experiments P-|]T5|, 



whose serious consideration for being associated with the primitive glueball is enhanced 
by the fact that its production is by mechanisms traditionally believed to be those that 
favour gluonic excitations. Specifically these include |16| 



1. Radiative J ftp decay: J/tp — > 7 + G [|17]] 

2. Collisions in the central region away from quark beams and target: pp —>■ Pf{G)p s 



12, 15 . 



Proton-antiproton annihilation where the destruction of quarks creates opportunity 



for gluons to be manifested. This is the Crystal Barrel |§-[|TTJ and E760 [13, [14| 
production mechanism in which detailed decay systematics of /o(1500) have been 
studied. 

4. Tantalising further hints come from the claimed sighting fl8| of the /o(1500) in 
decays of the hybrid meson candidate 7r(1800) — > 7r/ (1500) — > nr/r]. 

The signals appear to be prominent in decay channels such as 7777 and 777/ that are 
traditionally regarded as glueball signatures. However, such experiments are not totally 
novel and some time ago one of us (FEC) addressed the question of why glueballs had 
remained hidden during 25 years of establishing the Particle Data Group list p] of QQ 



states. This was suggested [16| to be due to the experimental concentration on a restricted 
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class of production mechanisms and on final states with charged pions and kaons. The 
more recent emphasis on neutral final states (involving 7r°, n, rf) was inspired by the 
possibiblity that 77 and n' are strongly coupled to glue. This dedicated study of neutrals 
was a new direction pioneered by the GAMS Collaboration at CERN announcing new 



states decaying to 7777 and r/77' [19 



The Crystal Barrel collaboration at LEAR has made intensive study of pp annihilation 
into neutral final states involving tt°, rj and n' . They find a clear signal for / (1500) in 
7r°7r°, 7777 and 7777'. Our present work extends and generalises the work of ref. JTBJ in the 
light of these new data from LEAR. Our purpose is to examine the data on the /o(1500), 
compare with predictions for glueballs and identify the seminal experiments now needed 
to confirm that gluonic degrees of freedom are being manifested in this region. A summary 
of this work has already been published elsewhere . 



The structure of the paper is as follows. We first review the experimental data on scalar 
mesons with special emphasis on states seen in the Crystal Barrel detector at LEAR. We 
then derive from SU(3)/ the branching ratios for QQ decays into two pseudoscalars, 
show that this successfully describes the known decay rates in the 2 ++ QQ nonet, and 
then compare our predictions to the observed decay modes of /o(1500). Previous bubble 
chamber experiments have not observed a KK signal in the 1500 MeV mass region which, 
if confirmed, would imply a set of branching ratios that are unnatural for a state belonging 
to a quarkonium nonet. If a significant signal were to be observed, it would be possible 
to find a quarkonium mixing angle that reproduces the observed final state abundances; 
however the systematics would then imply that / (1500) is dominantly nn = (uu+dd) / 
This would have two immediate consequences: 

1. This would leave the /o(1370) state, which is also seen in 7979 annihilation with decay 
branching ratios and total width consistent with an nn structure |^, 22|| , isolated. 



2. With either the /o(1370) or /o(1500) assigned as the nn member, the orthogonal 
quarkonium in the nonet would have to be dominantly ss, hence probably heavier 
than / (1500) and decaying strongly into KK. Identification of this state is now 
imperative in order to complete the multiplet and discriminate among hypotheses. 

We then show that that the decay rates of /o(1500) are compatible with a glueball state 
whose mass lies between the nn and ss scalar quarkonium states, and whose nearby 
presence disturbs the glueball decays in a characteristic flavour dependent manner. In 
the climax of the paper we show that dynamics inspired by lattice QCD may be consistent 
with the data and we consider the implications for glueballs mixing with quarkonia in the 
1500 MeV range. We finally show that a reasonable nonet can be constructed with the 
remaining scalar mesons. 

2 Scalar Mesons in the Crystal Barrel 

The lowest lying ++ mesons, namely the isospin 1 = / (980) and the 1=1 a (980), 



have been assumed to be KK molecules [23, 24]. This belief is motivated by their strong 
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couplings to KK - in spite of their masses being at the KK threshold - and their small 
77 partial widths. For / (980), T 77 = 0.56 ± 0.11 keV 0. For a (980), one finds with 



a LEAR measurement of the relative branching fraction for ao decay to KK and 7771" [25[] 
the partial width T 77 = 0.33 ± 0.13 keV 0. Thus the 77 partial widths appear to be 
nearly equal, close to predictions for KK molecules (0.6 keV) and much smaller than for 
QQ states [ZB|. 

The nature of these states is likely to be illuminated soon at DA0NE |24| . If they are 
not simply QQ then the ++ QQ mesons need to be identified. A new J PC (I G ) = ++ (l _ ) 
meson, a (1450) — > nn, has been reported by the Crystal Barrel collaboration at LEAR 



2§. This state, with a mass of 1450 ± 40 MeV and a width of 270 ± 40 MeV, appears, 
together with a (980), in the nn S-wave in pp annihilation at rest into 777r°7r°. We shall 
show in section 8 that ao(1450) can be identified with the 1=1 member of the ground 
state ++ QQ nonet. In turn this and the if* (1430) set the natural energy scale for the 
scalar nonet. 

Recent data in pp annihilation at LEAR into t/ttV J25|, 3tt° g | and r/r/vr JTCj, |22 



require an / = scalar resonance in the range 1320 to 1400 MeV, decaying to 7r°7r° and rjrj. 
We shall use 1360 MeV as average mass but shall adopt the nomenclature of the Particle 
Data Group calling this state /o(1370). Its width varies between 200 and 700 MeV, 
depending on theoretical assumptions. For example, the 37r° data give T ~ 700 MeV |§, 



decreasing to 300 ± 80 MeV if the 700 MeV broad background structure |27|, [28[] centered 
at 1000 MeV in the -kit S-wave (called /o(1300) in the latest issue of the Particle Data 
Group) is introduced in the analysis. This is in good agreement with a coupled channel 
analysis of nn°n°, 3tt° and 777771"° which leads to a mass of 1390 ± 30 MeV and a width of 



380 ± 80 MeV p9[. The 77 width in this region is also consistent with it containing the 



5 P (QQ) state gg, |T 



There is rather general agreement that the ground state nn state is manifested here. 
The debate is one of detail on the relationship of the /o(1370) to the broad /o(1300), in 
particular as to whether these are two independent states or manifestations of a single 
state, and to what extent unitarity corrections are important p2] . This issue is peripheral 
to our main analysis which will rely only on the generally accepted association of the 
ground state nn as the seed for the phenomenology in the /o(1300 — 1370) region and 
that this is distinct from the /o(1500) state observed by the Crystal Barrel collaboration 
decaying to 7r°7r° || 0], nn [|1(], ^2j and nn' [HI]. This state was seen in pp annihilation into 



3tt°, 77777r° and 777/71"°, leading to six final state photons. The masses and widths observed 
in the three decay channels are consistent, giving the average: 

(m, T) = (1509 ± 10, 116 ± 17) MeV, 



while the coupled channel analysis |29| gives 1500 ± 10 MeV and the less precise but 
compatible width of 154 ± 30 MeV. It is possible that / (1500) has also been seen by 
MARKIII and DM2 in J/ip — > 7 + tttttttt, hitherto misidentified as 0~ + |17| and in 



J/tp — » 77T7r[[3t|. A resonance decaying to n n and r/r] was also reported by E760 in 
pp annihilation at higher energies at the Fermilab accumulator, with masses and widths 
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(1508, 103) MeV [0 and (1488, 148) MeV respectively. A spin-parity analysis is in 
progress 



The GAMS Collaboration at CERN |]19| reports a ++ 180 MeV broad resonance, 



/o(1590), decaying to rjT]', rjr] and 47r°, also observed in central production [|15| and by 



VES at Serpukhov in 777/: this might be the / (1500) state though the status of the 
7T7T branching ratio needs to be clarified.^] A strong coupling of / (1500) to pions would 
contradict it being primarily an ss state and, as we shall argue later, it is a candidate 
for a glueball mixed with the QQ nonet, where /o(1370) is dominantly nn, and a more 
massive ss remains to be identified. The /o(1500), clearly established in different decay 
channels and with detailed information on branching ratios to several channels, will form 
the fulcrum of our investigation. 

There are candidates for this ss state though their existence and/or ss assignment 
remain to be established. A ++ structure, /q(1525), with poorly known width (~ 90 
MeV) is observed to decay into K S K S in K~p interactions [07[| . This state requires 
confirmation from other experiments. The #(1690) (now known as /j(1720) J7|) is a 
candidate due to its affinity for KK and f]f] decays though its quantum numbers, ++ or 
2 ++ , are still controversial. 

An 1 = scalar with a width of 56 MeV, is observed at 1446 MeV by the WA91 collab- 
oration at CERN in pp central collisions [|l!| . It decays to four pions, dominantly through 
P°(tt + 7t~)p where the dipion is in a P-wave. This may be the same as /o(1500) produced 
in the second of the favoured glueball mechanisms (section 1), with its apparent small 
width being the result of interference between / (1370) and / (1500) In any event 
it does not detract from the qualitative observation that there are too many isoscalars 
observed in various production mechanisms for them all to be explained naturally within 
a QQ picture. The fact that there does not appear to be such copious activity in the 
1 = 1 and strange sectors adds weight to the suspicion that glueball excitation is affecting 
the 7 = spectrum. 

A substantial part of this paper will examine what the flavour content of two-body 
decays can reveal about the structure of the initial meson. Based on this analysis we 
shall argue that the Crystal Barrel /o(1500) has decay properties incompatible with a QQ 
state. In addition we shall show that a reasonable QQ nonet may be constructed with 
a (1450), / (1370), K£(U30), and an ss state above 1500 MeV, thus leaving / (980) and 
/o(1500) as exotic (not simply QQ) states. 



3 Quarkonium Decay Amplitudes 

Consider a quarkonium state 

\QQ) = cosa I rm) — sina| ss) (2) 

1 The 7T7T branching ratio is the largest deviation but is "not in contradiction" with the /o(1590) of 
GAMS and the / (1500) of Crystal Barrel being the same state |p^| 
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where 

nn = (uu + dd) / \/2. (3) 
The mixing angle a is related to the usual nonet mixing angle 9 |7J by the relation 

a = 54.7° + 9. (4) 

For 9 = the quarkonium state becomes pure SU(3) f octet, while for 9 = ±90° it becomes 
pure singlet. Ideal mixing occurs for 9 = 35.3° (-54.7°) for which the quarkonium state 
becomes pure ss (nn). 
In general we define 

n = cos0|nn) — sin0|ss) (5) 

and 

n' = sin0|nra) + cos</>|ss) (6) 

with = 54.7° + 9ps, where 9ps is the usual octet-singlet mixing angle in SU(3)/ basis 
where 

rj = cos(9 PS )\r]8) ~ sm(9 PS )\r}i), (7) 

V = sm(0ps)|77 8 )+cos(0ps)|?7i>- (8) 

The decay of quarkonium into a pair of mesons QQ — > M(Qq,j)M(qiQ) involves the 
creation of q^qi from the vacuum. If the ratio of the matrix elements for the creation of 
ss versus uu or dd is denoted by f\ 

_ Wis*) (Q) 

p ~ {o\v\ddy { ) 

then the decay amplitudes of an isoscalar ++ (or 2 ++ ) are proportional to 



(QQ\V\irn) = cosa 

(QQ\V\KK) = cosa(p - v^tana)/2 

(QQ\y\w) = cosa(l -pv / 2tana)/2 

(QQ\V\nn') = cosa(l + pV2tana)/2. 



(10) 



The corresponding decay amplitudes of the isovector are 

(QQ\V\KK) = p/2 
(QQ\V\nri) = l/y/2 
(QQ\V\nn') = l/y/2, 



2 we shall assume that (0\V\dd) = (0\V\uu) 
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and those for K* decay 



{QQ\V\K7t} = V3/2 
(QQ\V\K V ) = (V2p-1) 
(QQ\V\Krf) = {V2p + 1) 

(12) 



For clarity of presentation we have presented eqn. |10] , [TT1 and [12] in the approximation 
where rj = (nn — ss) / v2 and rf = (nn + ss) / y/2, i.e. for a pseudoscalar mixing angle 
8 PS ~ —10° (0 = 45°). This is a useful mnemonic; the full expressions for arbitrary 
rj, r)' mixing angles Ops are given in appendix A and are used in detailed comparisons 
throughout this paper. Exact SU(3)/ flavour symmetry corresponds to p = 1; empirically 
p > 0.8 for well established nonets such as 1 and 2 ++ |2l], |39| . 

The partial width into a particular meson pair MiMj may be written as 

Y l3 = Cl] \M l3 \ 2 x \F t] (q)\ 2 x p.s.(q) = 7 J x {F^q}] 2 x p.s.(q) (13) 

where p.s.(q) denotes the phase-space, Fij(q) are mo del- dependent form factors, Mij is 
the relevant amplitude (eqn. [10,11 or [12]) and Cjj is a weighting factor arising from the 



sum over the various charge combinations, namely 4 for KK, 3 for ttti, 2 for 777/ and 1 for 
T]f] for isoscalar decay (eqn. [10]), 4 for KK, 2 for 7rrj and 2 for ttt]' for isovector decay (eqn. 
[TTJ) and 2 for i^* decays (eqn. |T2|). The dependence of 7?. = Cij\Mij\ 2 upon the mixing 
angle a is shown in fig. |l]a for the isoscalar decay in the case of SU(3)/ symmetry, p = 1. 

We confront the above with data on the established 2 ++ nonet, determine a probable 
range of values for p, and then compare with / (1500) decays. 



4 Flavour Symmetry in Meson Decays 
4.1 QQ Decays: the 2 ++ nonet 

Exact SU(3)/ flavour symmetry requires the parameter p to be unity. To get a feel- 
ing for symmetry breaking in the QQ sector we have computed some of the expected 
branching ratios for the tensor mesons a2(1320), ^(1270) and 7^2(1430) decaying to two 
pseudoscalars and have compared them with data 0. The decay branching ratio B of 
a QQ state is proportional to the partial width (eqn. |13|). We use for the phase space 
factor p.s.(q) = q and for the for mf actor 

\F^q)\ 2 = 1 2£ xexp(-q 2 /S(3 2 ), (14) 

where I = 2 is the angular momentum in the final state with daughter momenta q and 
(3 ~ 0.4 GeV/c p9| . The ratios of the various partial widths are rather insensitive to 



choice among different successful descriptions of meson spectroscopy and dynamics. The 
detailed sensitivity to form factors is discussed in appendix B. 



S 



The pseudoscalar mixing angle 9ps has the empirical value -(17.3 ± 1.8)° ||D| (and 
hence <fi = 37.4°, somewhat removed from the "ideal" 45° used in section 3). We shall use 
the full expressions given in appendix A and the above value of 9ps in all phenomenology. 
Analogously, for the tensor mixing angle, 9 = 26° f7|. The predictions are fitted to the 
experimental values [|7| 



B(f 2 - VV)/B(h - tttt) = (5.3 ± 1.2) x 1(T 3 , 
B(f 3 -> KK)/B(f 2 -> tttt) = (5.4 ± 0.6) x 1(T 2 , 
B{a 2 -> V^)/B(a 2 -> = 2.95 ± 0.54, 
5(a 2 -> r/rf/Bfa = 0.116 ± 0.029. 



(15) 



Figure § shows the x 2 distribution as a function of p for various values of /3. The 
distribution does not change significantly for (3 >1 GeV/c. A good fit (with a x 2 confidence 
level of more than 5%) is obtained with (3 = 0.5 GeV/c (or larger) for which p = 0.96 
± 0.04. This result is consistent with 1^(1430) decays although the experimental errors 
are large: B(K* -> Kn)/B(K* -> Kr?) = (3.6t££) x 10 3 |]. We therefore conclude that 
flavour symmetry breaking effects cannot be large in this established QQ nonet. Similar 
conclusions follow for a wide range of QQ decays (at least in the (3 — > oo limit where form 
factors are ignored Thus it seems reasonable to expect that for QQ scalar decays 

also, p w 1 and /3 ~ 0.5 GeV/c. 

We have refrained from using the corresponding ratios for /2(1525) decays since the 
7T7T decay width is (i) poorly known and (ii) very sensitive to the precise tensor nonet 
mixing angle 9. For the ratio B(f^ — > r\r\) j B{f 2 — > KK) we find 0.07 for p = 1 and /3 = 
0.5 GeV/c, in good agreement with experiments (0.11 ± 0.04) |7j but at variance with the 
value (0.39 ± 0.05) advocated by the Particle Data Group which relies on one experiment 
only. 



4.2 The Decay Properties of / (1500) 

The branching ratios for /o(1500) production and decay are 



B(pp - f n°, f -> 
B(pp^ f 7t°J 
B(p P ^f 7T°J 



^0 On 
71 71 



Wl) 

vv r ) 



(8.1 ±2.8) x 10" 4 , 
(5.5 ±1.3) x 10" 4 , 
: (1.6 ±0.4) x 10 -4 . 



11, 22 



(16) 



where the errors do not reflect the statistical significance of the signals, but rather uncer- 
tainties in the various assumptions made in the fitting procedures. The decay branching 

q and the form factors 



ratios are given by eqn. 13 with p.s.(q) : 



exp(-g 2 /8/3 2 ), 



(17) 
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since t = for ++ decays to two pseudoscalars. The /o(1500) observed in nn' decay has 
a mass of 1545 ± 25 MeV and lies just above threshold jnj. We use for q the average 
decay momentum (194 MeV/c) derived from the damped Breit-Wigner function used in 
the analysis of ref. |]TT[. The uncertainty in the mass (± 25 MeV) is taken into account 
when computing the error on 7 2 . With f3 = 0.5 GeV/c (section 4.1) we find: 

Rl = fiMm^4 = 0.27 ±0.11, (18) 



7 2 (/o(1500) 



7T7T 



ft= ffiS^ =0.19±a08, (19) 
7 2 (/o(1500) -> tttt) 

where 7T7t includes 7r + 7r~. These results are in good agreement with the results of the 
coupled channel analysis |2{| . A signal for scalar decay to KK has not yet been observed 



in pp annihilation in the the 1500 MeV region. A bubble chamber experiment j41] reports 
B.R.(pp —> Xtc,X — > KK) < 3.4 x 10 -4 which interpreted directly as an intensity leads 
to a (90% C.L.) upper limit: 

7 2 (/o(1500)^A-7r) < 

* = fUtisoo)^™) °' L (20) 

However, interference effects among amplitudes could lead to an underestimate for this 
number. We shall consider the implications of the above R3 value but shall also give 
results allowing for larger values. 

We find that if in the decay of some state the ratios of partial widths (per charge 
combination and after phase space and form factor corrections) for nn/rnr and rjrj/KK 
are simultaneously both greater than unity, then this state cannot be a quarkonium decay 
unless ss production is enhanced (p > 1) (see fig. ^). The / (1500) data on rjrj/KK 
satisfy this, but the 7777/ 7T7t is inconclusive; at la the ratio per charge configuration gives 
0.81 ±0.33. If / (1590) and / (1500) are the same state, as discussed above (see footnote 
1), then if the tttt branching ratio is reduced towards the GAMS limit ||T5|, |I~9|j the value of 
Ri would rise such that it may be possible to confirm the /o(1500) as a glueball by this 
test alone. We cannot overemphasise the importance of a mutually consistent analysis of 
the data on these experiments, in particular for clarifying the magnitude of the nn/nn 
ratio. If the ratio rises, as for GAMS, it would immediately point towards a glueball; if the 
ratio remains as in eqn.|l8] then the arguments are less direct but there still appears not 
to be a consistent QQ solution to the flavour dependence of the ratios of partial widths 
and the magnitudes of the total widths for the / (1500) — / (1370) system. 

Figure |I] shows the invariant couplings j 2 as a function of a for p — 1 (a) and p = 0.75 



and 1.25 (b), for a pseudoscalar mixing angle 6 = -17.3° |j4"0"f . The effects of SU(3)/ 
breaking, p < 1, in the region where nn dominates the Fock state ( < a < 30°) are 
interesting (fig. |l|b). We see that the branching ratios for 77 or v( are little affected (essen- 
tially because they are produced via the nn component in the 77 which is p independent) 
whereas KK depends on p (due to ss creation triggering KK production from an nn 
initial state). Thus we can suppress KK by letting p — ► without affecting the nnj-n-n 
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ratio substantially. In this case the measured values for nn and nn' and the upper limit 
for KK suggest that a ~ 0, hence /o(1500) is a pure nn meson. However this is still 
unsatisfactory as the required value of p implies a dramatic suppression of ss creation to 
a degree not seen elsewhere in hadron decays. 

Figure [| shows the dependence of p tana (eqn. |10|) on Ri and R 2 - From the experi- 
mental values (eqn. |T8]and |TU|) we find consistency for p tan(a) ~ -0.1. Hence either p is 
small or a is small. The former leads to unacceptable violation of SU(3)/ and the latter 
predicts the ratio j 2 (KK) /^(tctt) to be 1/3. 

Form factors of the type in eqn. [17] tend to kill transitions at large q, for example 
F w (g^)/F 7r7r (g 7r ) > 1 (i.e. opposite to naive phase space which grows with q ). Without a 
form factor, e.g. = 1, we obtain R x = 0.31 ±.13, R 2 = 0.25 ±0.11, which excludes 

a common range of p tan(a) (see fig. and R3 < 0.12. The form factors used in ref. 
|1| have a node at q ~ 0.9 GeV/c and hence lead to an even stronger suppression of the 



observed Tin intensity which dramatically reduces R3, in contradiction with the expected 
1/3 for an nn state (see also appendix B). 

In fig. [5] we plot the allowed regions of p vs. a. The grey area shows the common 
values of p and a which satisfy the Crystal Barrel data each at the 90 % C.L. while the 
black area shows the restricted range allowed by KK (eqn. 20). 

If one wishes to force /o(1500) into a QQ nonet, then independent of form factors and 
SU(3)/ breaking one is forced to a — > 0, whereby /o(1500) has strong nn content. This 
remains true even were the KK branching ratio, currently being remeasured at LEAR, 



significantly greater than the R% value of eqn. 20: the magnitude of the KK/tctt ratio 



is controlled more by flavour symmetry breaking than by the magnitude of the nn - ss 
mixing angle in the a — > region. This immediately implies that the orthogonal isoscalar 
member will be dominantly ss, with mass above 1500 MeV and prominent in KK. The 
/o(1525) if confirmed or the "0" (/j(1720)), if ++ , could be this state. However, the 
/o(1370), seen in both 7nr,rjn is then left in isolation. 

In the next sections we confront the data on G = /o (1500) with the extreme hypothesis 
that it is dominantly a glueball. Its production in the canonical glue enhanced environ- 
ments of J/V> -> 7(G -> 4vr) 0, pp -> p(G)p H, and pp annihilation g, g jTO], 0, 0, [l| 



is consistent with this hypothesis, and recent lattice QCD studies [||, [| predict that a 
scalar glueball exists in this region of mass. It thus scores well on two of the three glueball 
figures of merit [[16]. 

We shall now consider the dynamics and phenomenology of glueball decays. 



5 Primitive Glueball Decays 

The decays of cc, in particular Xo,2, provide a direct window on G dynamics in the 
++ ,2 ++ channels insofar as the hadronic decays are triggered by cc — > gg — > QQQQ 
(fig. |^a). It is necessary to keep in mind that these are in a different kinematic region 
to that appropriate to our main analysis but, nonetheless, they offer some insights into 
the gluon dynamics. Mixing between hard gluons and ++ , 2 ++ QQ states (fig. ||c) is 



11 



improbable at these energies as the latter 1-1.5 GeV states will be far off their mass- 
shell. Furthermore,the narrow widths of Xo,2 are consistent with the hypothesis that the 
3.5 GeV region is remote from the prominent + ,2 + glueballs, G. Thus we expect that 
the dominant decay dynamics is triggered by hard gluons directly fragmenting into two 
independent QQ pairs (fig. |6]a) or showering into lower energy gluons (fig. |^b). We 
consider the former case now; mixing with QQ (fig. 6c) and G — > GG (fig. 6b) will be 
discussed in section 6. 

G -> QQQQ 



This was discussed in ref. [I6"| and the relative amplitudes for the process shown in fig. 
|]a read 

(G\V\irir) = 1 
(G\V\KK) = R 
(G\V\ VV ) = (l + R 2 )/2 
(G\V\rrf) = (l-R 2 )/2, 

(21) 

with generalizations for arbitrary pseudoscalar mixing angles given in appendix A and 
where R = (g\V\ss)/{g\V\dd). SU(3)/ symmetry corresponds to R 2 = 1. In this case the 
relative branching ratios (after weighting by the number of charge combinations) for the 
decays Xo,2 - * Wi VV'i would be in the relative ratios 3 : 1:0:4. Data for xo 
are in accord with this where the branching ratios are (in parts per mil) [^]: 



\b{k+k-) 



3.1 ±0.6 
3.7 ±1.1 

3.5 ± 1.2 
2.5 ± 1.1. 



(22) 



No signal has been reported for r\vj '. Flavour symmetry is manifested in the decays of 
X2 also: 



B 7T U 7T 



-A 



1.1 ±0.3 
0.95 ±0.50 



-B(K + K-) 
B(rm) 



0.75 ±0.55 
0.8 ±0.5, 



(23) 
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again in parts per mil. The channel nn' has not been observed either. These results 
are natural as they involve hard gluons away from the kinematic region where G bound 
states dominate the dynamics. If glueballs occur at lower energies and mix with nearby 
QQ states, this will in general lead to a distortion of the branching ratios from the 
"ideal" equal weighting values above (a detailed discussion of this follows in section 6.1), 
and also in causing significant mixing between nn and ss in the quarkonium eigenstates. 
Conversely, "ideal" nonets, where the quarkonium eigenstates are nn and ss, are expected 
to signal those J PC channels where the masses of the prominent glueballs are remote from 
those of the quarkonia. 

An example of this is the 2 ++ sector where the quarkonium members are "ideal" . Data 
on glue in the 2 ++ channel, and potential mixing of glue with nn/ss, may be probed by 
J/ij) — > 7 + /2( 1270) //2(1525) which measures the gg — > nn/ss amplitude (fig. ||a) insofar 
as J/ip — > 7 + gg mediates these channels. The branching ratios in parts per mil are 

-B{J/ijj -»• 7/2(1270)) = 0.69 ±0.07 
^(7/^^7/2(1525)) = 0.63 ±0.1. 

(24) 

Here again there is no sign of significant symmetry breaking. Furthermore we note 
the ideal nn and ss nature of the 2 ++ , manifested both by the masses and the flavour 
dependence of the branching ratios, which suggests that G mixing is nugatory in this 
channel. These data collectively suggest that prominent 2 ++ glueballs are not in the 
1.2 — 1.6 GeV region which in turn is consistent with lattice calculations where the mass 
of the 2 ++ primitive glueball is predicted to be larger than 2 GeV. The sighting of a 
2 ++ state in the glueball favoured central production, decaying into nn with no significant 
7T7rf42|| could be the first evidence for this state. In view of our earlier remarks on the nn/mr 
and nn/ KK ratios being a potentially direct signature for a glueball, we recommend that 
a detailed search now be made for this state in nn and KK (and ittt) channels in central 
production. 

The phenomenology of the J PC = ++ sector in the 1.2 - 1.6 GeV region is rather 
different to this: the /o(1500) - /o(1370) system cannot be described within a QQ nonet, 
nor do the decay branching ratios of the / (1500) respect the flavour blindness of glue, 
(eqn. 19 and 20). We shall now begin to focus on this problem. 

It was shown earlier jnj that violation of flavour symmetry (R 2 7^ 1) leads to smaller 
KK/nn and a finite nn', at least if graph ||a dominates glueball decay. This follows 
immediately from eqn. [H] or from the generalized formulae given in appendix A. The 
contributions of graph ||a are shown in fig. |7|a as a function of R 2 . Graph |^a becomes 
compatible with the Crystal Barrel data and the small KK ratio if \R\ ~ 0.3 0] - a rather 
strong violation of symmetry which might be suggestive of significant mixing between G 
and flavoured states. If R3 were as large as | (the value for an nn state) an attempt to 
interpret as gluonium would still require \R\ to be as small as 0.5. We now show that 
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mixing of G and QQ is to be expected if the strong coupling picture of QCD, as in the 
lattice, is a guide to their dynamics. 

6 QQ and Glueball Decays in Strong Coupling QCD 

In the strong coupling {g — > oo) lattice formulation of QCD, hadrons consist of quarks 
and flux links, or flux tubes, on the lattice. "Primitive" QQ mesons consist of a quark 
and antiquark connected by a tube of coloured flux whereas primitive glueballs consist of 
a loop of flux (fig. ||]a,b) [[|3|]. For finite g these eigenstates remain a complete basis set 



for QCD but are perturbed by two types of interaction [44]: 



1. V\ which creates a Q and a Q at neighbouring lattice sites, together with an ele- 
mentary flux-tube connecting them, as illustrated in fig. |8]c, 

2. Vi which creates or destroys a unit of flux around any plaquette (where a plaquette 
is an elementary square with links on its edges), illustrated in fig. |8|d. 

The perturbation V\ in leading order causes decays of QQ (fig. ||s) and also induces 
mixing between the "primitive" glueball (Go) and QQ (fig. |8)f). It is perturbation V2 in 
leading order that causes glueball decays and leads to a final state consisting of G0G0 (fig. 
Da;); decays into QQ pairs occur at higher order, by application of the perturbation V\ 
twice. This latter sequence effectively causes G mixing with QQ followed by its decay. 
Application of V 2 leads to a Q 2 Q 2 intermediate state which then turns into colour singlet 
mesons by quark rearrangement (fig. |B|a); application of V2 would lead to direct coupling 
to glue in r], rf or V2 x V 2 to their QQ content (fig. ||b). 

The absolute magnitudes of these various contributions require commitment to a de- 
tailed dynamics and are beyond the scope of this first survey. We concentrate here on their 
relative contributions to the various two body pseudoscalar meson final states available 
to ++ meson decays. For QQ — > QqqQ decays induced by V\, the relative branching ratios 
are given in eqn. |10] where one identifies 

{QssQWQQ) 
P (QddQ\Vi\QQ)' 1 ' 

The magnitude of p and its dependence on J PC is a challenge for the lattice. We turn 
now to consider the effect of V± on the initial "primitive" glueball Go- Here too we allow 
for possible flavour dependence and define 

~ (ddWGo)' [ ] 

The lattice may eventually guide us on this magnitude and also on the ratio R 2 / p. In 
the absence of this information we shall leave R as free parameter and set p — 1. 
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6.1 Glueball-QQ mixing at 0{V X ) 

In this first orientation we shall consider mixing between G (the primitive glueball state) 
and the quarkonia, nn and ss, at leading order in V\ but will ignore that between the two 
different quarkonia which is assumed to be higher order perturbation. 
The mixed glueball state is then 

a =\G ) + {nn\Vi\G ) |ss) (ss\Vi\G ) , . 
which may be written as 

G = \G ) + ^ lVllG °J A V2\nn) + coR*\ss)} (28) 

where 

w = ~ ^l, (2Q) 

is the ratio of the energy denominators for the nn and ss intermediate states in old 
fashioned perturbation theory (fig. ||d). 

Denoting the dimensionless mixing parameter by 

_ (ddlV^Gp) 

the three eigenstates become, to leading order in the perturbation, 



iV G |G> = | Go) + ^{v^lnn) + rf|ss>} = |G > + v^f |QQ) 

iV s |* s ) = |ss)-^|G ) 
7V„|^ n ) = |nn)-eV2|G ) 



with the normalizations 



N G = Jl + ?{2 + u*R% 



N n = Jl + 2?. 



(31) 



(32) 



Recalling our definition of quarkonium mixing 

\QQ) = cosa|nn) — sina| ss) (33) 
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we see that Go has mixed with an effective quarkonium of mixing angle a where -y/2tana = 
—ujR 2 (eqn . [K|). For example, if uoR 2 = 1, the SU(3)/ flavour symmetry maps a glueball 
onto quarkonium where tana = — l/v2 hence 9 = —90°, leading to the familiar flavour 
singlet 

\QQ) = \ U u + dd + ss) /V3. (34) 

When the glueball is far removed in mass from the QQ, uj — > 1 and flavour symmetry 
ensues; the Xo,2 decay and the 2 ++ analysis of sections 4 and 5 are examples of this "ideal" 
situation. However, when u ^ 1, as will tend to be the case when Go is in the vicinity of 
the primitive QQ nonet (the ++ case of interest here), significant distortion from naive 
flavour singlet can arise. 

If the Go component contributed negligibly to the decays, the expectations would be 
that there is (eqn. |31]) (i) a state ^ n — > tt7c,i]1], KK which is compatible with / (1370); 
(ii) a state ^> s — > KK,7]r],7]ri', but not tttt, to be established; (iii) the 1500 MeV state G 
for which the decay amplitudes relative to nil are (replacing \/2 tana by —ujR 2 in eqn. 
0), 



(G\V\nn) = 1 
(G\V\KK) = {p + uR 2 )/2 
(G\V\ VV ) = (1 + uj P R 2 )/2 
(G\V\ VV ') = (l-u P R 2 )/2, 

(35) 

with generalisation given in appendix A. The invariant decay couplings 7^ are shown in 
fig. 0c as a function of uR 2 for p = 1. Thus, for example, SU(3)/ may be exact for the 
glue-quark coupling (R = 1) but mass breaking effects (Am = m s — 7^ 0) can cause 
dramatic effects if Eg — E nn or Eg — E sS is accidentally small, such that uj — > or 00 
respectively. We now consider an explicit mixing scheme motivated by three mutually 
consistent phenomenological inputs: 

1. The suppression of KK in the /o(1500) decays suggests a destructive interference 
between nn and ss such that ujR 2 < (see fig. ^c). This arises naturally if the 
primitive glueball mass is between those of nn and the primitive ss. As the mass 
of Go — > m nn or m sS , the KK remains suppressed though non-zero; thus eventual 
quantification of the KK signal will be important. 

2. Lattice QCD suggests that the "primitive" scalar glueball Go lies at or above 1500 
MeV, hence above the 1 = 1 QQ state a (1450) and the (presumed) associated nn 
/o(1370). Hence E Go — E nn > in the numerator of to. 

3. The Am = m sS — m nn m 200 — 300 MeV suggests that the primitive ss state is in 
the region 1600-1700 MeV. This is consistent with the requirement from (1) and (2) 
that m nn < m Go < m sS . 
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Higher order perturbation effects will be required for a complete treatment, in partic- 
ular including mixing between nn and ss, but that goes beyond this first orientation and 
will require more data to constrain the analysis. We shall present a posteriori evidence 
supporting this leading order approximation. 

Tests of this scenario and its further development will follow as the predicted states 
are isolated and the flavour dependence of their branching ratios is measured. In order to 
compute the decay branching ratios of the physical (mixed) states, we need to incorporate 
the contributions from the primitive glueball components, Gq. We consider this now. 



6.2 G -► G G at 0(V 2 ) 

Here the glueball decays directly into pairs of glueballs or mesons whose Fock states have 
strong overlap with gg (fig. ^b). This topology will not feed final states such as tttt nor 
KK since gluons are isoscalar. To the extent that there is significant G coupling to 77, rf 
or to the tttt S-wave, (mr) s , (e.g. ip' — ► ipr) and iP(tttt) s each have large intrinsic couplings 
notwithstanding the fact that they are superficially OZI violating) one may anticipate 
7777,777/, and (7T7r) s (7r7r) s in the decays of scalar glueballs. Analogously for CT + glueballs 
one may anticipate r](7nr) s or i]'(tttt) s decays. 

The manifestation of this mechanism in final states involving the 77 or 77' mesons 
depends on the unknown overlaps such as (gg\V\qq) in the pseudoscalars. We consider 
various possibilities from the literature without prejudice at this stage. 

In the limit m u ^ — > chiral symmetry suggests that the direct coupling of glue to 
the 77 or 77' occurs dominantly through their ss content, thereby favouring the 77'. This 
argument has been applied to the 7/(1460) in refs. p5| pj , [47jj : 



(gg\V\rf) _ (ss\ V ') _ cot{4>) + v / 2A_ 



(gg\V\r)) (ss\r]) A v / 2cot(0) 
where 



(gg(0-)\V\ss) 

in the chiral limit, for which the ratio in eqn. ^ is ~ —4/3. 



(37) 



The ratio eqn. |36] depends sensitively on the pseudoscalar mixing angle and on a small 
breaking of chiral symmetry but remains negative in the range —0.9 < A < 0.5. Thus we 
anticipate 

_ ( VV >\V\G ) (gg\V\rf) _4 
r °" (r777|y|Go) (gg\V\ V ) 3" 1 ' 

There is some ambiguity as to how this is to be applied quantitatively since my 7^ 
m v and the wavefunctions at the origin 7/y(0) and 7/^(0) are in general different. An 
alternative measure [|5^] may be the ratio r(-^ — > , yr] f )/r(ip — > 777) = 5.0 ± 0.6. Dividing 
out phase space factors ~ pi we obtain the ratio of matrix elements 



r (J/i>) = ±(2.48 ± 0.15). (39) 
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The solution with the negative sign is compatible with a small breaking of chiral symmetry 
(A = 0.18). This gives similar results to arguments based on the gluon anomaly in the 



pseudoscalar channel (see eqn. (60) in ref. ||8| and also ref.[j49 



This enhanced gluonic production of 7777, 777/ does not appear to be dramatic in the 
Xo,2 decays as the 7777 / 7T7t ratio appears to be "canonical" in the sense of section 5. This 
may be because (gg\V\rj) and {gg\V\r]') are hidden in the large errors on the present 
data (in which case isolation of \ — > rj'rj' at a Tau Charm Factory would be especially 
interesting). Alternatively it may be (<7<7|V|(7T7r) s ) that is important and hence x ~^ 47T 
is the signal. Indeed this channel is the biggest hadronic branching ratio for both xo and 
X2- It would be interesting to compare the tttt spectra in these final states with those in 
Vc Vjirn and 777r7r which are dominant modes in the r\ c decays and may also be signals 
for this dynamics. High statistics from a Tau-Charm Factory may eventually answer 
this question. The relative coupling strengths of rjim and rj'im in decays of the glueball 
candidate 7/(1420) are also relevant here. 

Note that for the decays of J PC = ++ , 2 ++ . . . states, one will have a sharp test for a 
glueball if non-perturbative effects favour the direct 7777, 7777' decay path over mixing with 
QQ systems. A state decaying to 7777, 7777' and/or 7/7/ but not tttt nor KK cannot be simply 
QQ since p and a cannot sensibly suppress mr and KK simultaneously - (see e.g. fig. 
and D. GAMS has claimed states at 1590 MeV (7777,7777') P, 1740 MeV (ijrf) H 
and 1910 MeV (7777') ]5T] with no strong signal seen in tctt nor KK. If the existence of 
any of these enigmatic states X is seen in other experiments such as central production 
pp — > p{X)p or pp annihilation, where X — > rjr), rjr)' with no KK signal, this will be strong 
evidence for the presence of Go in their wavefunctions. The possible hints of /j(2100) in 77 
channels |T^] and of ^(2170) in 7777 [|2|, if confirmed, will put a high premium on searching 
for or limiting the KK branching ratios for these states. A particular realisation of these 
generalities is the model of ref. ||52|| . Similar remarks apply to the decay of pseudoscalars 
— ► 77(7T7r) s in contrast to K(Kn) s . This may be the case for the 7/(1440) and for the rj c , 
as discussed above. 



7 Application to Scalar Mesons around 1.5 GeV 

7.1 Decays of / (1500) 

We shall now combine these ideas with the other result of section 6, namely that at 0(Vi) 
the Go mixes with QQ with amplitude £ and that the resulting QQ components decay as 
in section 3. 



For simplicity of analysis we shall set p = 1. From eqn. |3T] and [35| we obtain for 
G = / o (1500): 

(KK\Vx\G) = 1+ ^ R \ nn\V l \G). (40) 

Eventual quantification of R 3 may be translated into a value of uR 2 (see fig. §). We shall 
scale all decay amplitudes relative to that for (7T7r|V|G) and see what this implies for the 
G decay amplitudes. Thus 
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from Ri (eqn. [18]) and 



from i?2 (eqn. |19j) , from which we predict the ratio 



(VV'\V\G ) = 2r 2 - sin20(l - uR 2 ) 
r ° (w\V\G ) 2n - cos0(l - cui? 2 ) - 1 - wi? 2 ' 1 j 



The ratio rg is plotted in fig. [TO as a function of ujR 2 . Note that this ratio is rather 



sensitive to the precise value of the mixing angle 0. The solution with the + sign in 
eqn. ^1] and the — sign in eqn. |42| (which we refer to as the "+ -" solution) agrees very 



well with ujR 2 in the range predicted by the current KK suppression (fig. |9|) and with 
radiative J /if) decay (eqn. |39[). Figure |10| also suggests another possible solution ("- -") 
compatible with the if) — > 777/77' ratio. 

In the particular limit ujR 2 — — 1, the 7777 decay mode would be driven dominantly by 
Go decay. The smaller rate for 7777' decay observed by Crystal Barrel would then be due to 
destructive interference between Go decay and the admixture of quarkonium in the wave 
function common to the "+ -" and "- -" solutions. The amplitudes for Go decay (fig. []b) 
are given in appendix A and the invariant couplings shown in fig. [7]b as a function of A. 

7.2 / (1370) 

The decay amplitudes for \I/ n — > tttt and KK will be those of an nn state such that, 

7 2 (^„ -> KK) 1 



7 2 (^„ -> tttt) 3 



(44) 



Assuming that the lower state \& n is /o(1370) the decay to r/r]' is kinematically forbidden 
and so the 7777 decay will be the only one immediately sensitive to the predicted Go 
component in the /o(1370) Fock state. We find from the decay branching ratios measured 
by Crystal Barrel @, |22| Q 



B(pp - /o(1370)7r°, / - vrV) < (2.6 ± 0.4) x 10" 3 , 
B(pp -> /o(1370)7T°, /o -> rjrj) = (3.5 ± 0.7) x 10" 4 , 

(45) 



3 The < sign reflects the fact that the measured branching ratio for /o(1370) decay to tttt also includes 
some contribution from /q( 
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after phase space and form factor corrections: 



7 2 (^ n -> rm) 



or 



(^ n |y|7T7r) 



> 0.07 



> 0.46. 



On the other hand, eqn. 



IT] predicts 



cos 



V 7 ^ 



(G |V|w> 



(46) 
(47) 

(48) 



which reduces to cos 2 <; 
into eqn. f|I] one finds with eqn. [47 



(^ n |y|7T7r) "" v ~^{nn\V\^Y 

0.63 for a pure rtn state. Solving for (Go\V\rjrj) and introducing 



cos 



0.46 



2n - (1 + cui? 2 ) - cos 2 0(l - u;i? 2 ) 



(49) 



The result eqn.[|7] then implies that |£| must be small and prefers the +- solution rather 
than the - - solution. One finds the 90% CL. upper limit 



If I < 0.47 



(50) 



which a posteriori justifies the first order perturbation used in the derivation of eqns. ^7 
and 0. 

Further data analysis is now needed to quantify the experimental ratio eqn. |46| and 
compare it in detail with the above. In any event, the branching ratio of this /o(1370) state 
is consistent with nn dominance and hence further isolates the /o(1500) as an exceptional 
state. 



7.3 The / (1370) - / (1500) system 

It should by now be clear that it is the combination of the two siblings, /o(1370) and 
/o(1500), rather than either one on its own that reveals the need for degrees of freedom 
beyond qq. We now illustrate how the branching ratios and widths of the pair manifest 
this quantitatively. 

In the quark model of ref. B§[ the widths of 3 Pq are qualitatively ordered as T(nn) > 



T(ss) > T(a ) > T(K*). Empirically T(a ) = 270 ± 40 MeV, F(K*) = 287 ± 23 MeV 
which supports this pair to be members of the nonet and leads one to expect for their 



partners that T(nn) ~ 700 MeV and T(ss) ~ 500 MeV. In the flux tube model of ref.p4 
(see also appendix B) after normalising to the known widths of the 2 ++ nonet one finds 
typically T(Kq) > 200 MeV and T(a ) > 300MeV in accord with data, and predict 
r(/,5 in ) > 500MeV. That the /™ width will be very broad is a rather general conclusion 
of all standard quark models (see also ref. ||53|| ); unitarisation effects do not alter this 
conclusion 
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The /o(1500) width of 116 ± 17 MeV is clearly out of line with this, being even smaller 
than the K* and ao widths. The total width of /o(1370) is not yet well determined, 200- 
700 MeV being possible || ^] depending on the theoretical model used in the analysis. 



These suppressions of widths are natural in the G — QQ mixing scheme as the presence 
of the Go component dilutes the effect of the leading nn component: 

r(* n ) = r(rm)/(l + 2£ 2 ) > ^r(nn). (51) 

The /(1500) by contrast is, in our hypothesis, a glueball in leading order and with nn, ss 
components at 0(£) in perturbation. The decay amplitudes for /(1500) are all at 0(£), as 
shown by our analysis above. 

Quantitative measures arise if we concentrate on the decays into pseudoscalar pairs. 
The measurements from Crystal Barrel 22] give the ratios of branching ratios for 
G = / (1500) as 

^ = 0.23, ^ = 0.07 (52) 

TTTT TTTT 

which implies that B(G — > tttt) = 0.7F 2 G (where is the fraction of of two-body decays). 
With Tq = 116 MeV this implies that Y{G — > tttt) = 28.2i<2 MeV per charge mode and 
hence, after dividing out phase space and form factors (738 MeV/c and 0.755 respectively) 
we have the reduced dimensionless measure 

f(G -> tttt) = i 7 2 (G -> tttt) = 0.05F 2 G . (53) 

We now perform the same manipulations for the / (1370). The tttt decay appears to 
be the dominant two-pseudoscalar decay mode (eqn. |45]) ; any error on neglecting nn in 
the analysis is likely to be masked anyway by the uncertainty on the total width which 
is not known to better than a factor of 3.5 (spanning 200 - 700 MeV). Dividing out the 
phase space (671 MeV/c) and form factors (0.798) as before, we form the reduced measure 
per charge combination 

f(/(1370) -> tttt) = \-i 2 {G -> tttt) = (0.13 - 0.44) F/. (54) 
Hence the ratio of measures for the two states is 



no ^ tttt) 



G 
2 



T(f ^ tttt) F, 



< OA-*?. (55) 



f 
2 



Assuming that R = 1 but allowing u and £ to be free, we expect this ratio to be given 

by 

f (g -> tttt) 2e(i + 2e) 

f(f^TTTT) l + e(2 + 00 2 ) { ) 

If Fl is not small, then independly of F 2 G : 
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If I < 0.52, (57) 



which is consistent with our earlier, independent, estimate in eqn. [50[ The sum of the 
partial widths of the two states for the tttc channels is 

r(*„) +T(G) ~ T(nn). (58) 



Ref. |54] reports a strong ~ 300 MeV 47r signal in the 1400 MeV region. It is unlikely 
that this signal is /o(1370) since the inelasticity in the 7T7T S-wave would be very large (~ 
80 %). However, given the uncertain dynamics in the irn sector one must allow for this 
possibility in which case our result eqn. [57| would break down. On the other hand, a Air 
contribution to /o(1500) decay would decrease the upper limit for |£|. A more detailed 
analysis is now warranted to verify if this further qualitative indication is supported and 
to quantify the resonant contributions F^/F^ and their effect on the analysis above. 

7.4 ^("1600") 

The mass of the state, the value of ujR 2 and of m(Go) are all related in our scheme and 
currently we can say no more than that the 1520 MeV to 1850 MeV range is possible for 
the mass of the m(\I/ s ). A small value of |£| would then require that \fr s decays essentially 
like an ss state, with couplings 

7 2 (7T7r : KK : rfq : rpf) = : 4 : - : 2. (59) 

The decay of / (1500) to KK will set the scale. The amount of G mixing depends 
rather sensitively on the value of u. However, a general result is that there will be 
constructive interference between Go and ss for the r]T]' channel in either the "+-" or "- 
-" solutions 

7 2 (fr, -> rpf) j 2 (ss -» ryrf) 1 

^(^s^KK) >y*{8B^KK) 2' 1 ' 

If we take the widths of if* (1430) and ao(1450) as a guide to normalize the width of 
the nonet in the Godfrey-Isgur model [j39~| , we would anticipate T(ss) ~ 500 MeV at a 
mass of 1600 MeV. The \l/ s width will be suppressed relative to that of a pure ss due to 
the glueball component (reflected in the normalization, eqn. |31|) but the actual branching 



ratios may be sensitively dependent on dynamics. In the model of ref. jPfl the KK and 
i]!] widths are suppressed if m(^ s ) — > lSOOMeV", see appendix B. The importance of the 
f]i]' channel appears to be a solid prediction. 

If the /j(1710) is confirmed to have a J = componenet in KK but not in tttt, this 
could be a viable candidate for a G - ss mixture, completing the scalar meson system 
built on the glueball and the quarkonium nonet. 
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Table 1: The scalar mesons, their observed decay modes and their suggested assignment. 



Isospin 


State 


Decays 


Nature 


1 


a (980) 


rjn, KK 


KK molecule 





/o(980) 


7T7T, KK 


KK molecule 


1 


a (1450) 


1]7l 


QQ 0++ nonet 





/o(1370) 




Q~Q 0++ nonet 




[/o("1600")] 


[KK, 7]r), rjr)'] 


QQ ++ nonet 


1/2 


fT *(1430) 


Ktt 


QQ ++ nonet 





/o(1500) 


7T7T, T]T], T]T]', 47T° 


glueball 



8 The Scalar QQ Mesons 

Based on the analysis of the previous sections we shall assume that /o(1500) is mainly 
glue and shall examine whether a reasonable scalar QQ nonet can be constructed with 
the other scalar mesons, neglecting first the small glue admixture in the two mainly QQ 
isoscalars. 

As we have seen, there are too many isoscalar ++ mesons to fit in the QQ ground 
state nonet. A possible classification of the scalar mesons is shown in table [l]. The a (980) 
and /o(980) are interpreted as KK molecules, a hypothesis that may be tested at DA$NE 
IJ. The a (1450), / (1370), /q(1600) and fsT *(1430) are the members of the ground state 
QQ scalar nonet. Note that the masses of the strange and 1=1 members are similar; 
this is also the case for some other nonets, in particular 4 ++ 0. 

The Crystal Barrel and Obelix Collaborations at LEAR also report the observation of 
a ++ state decaying to p + p _ and p°p° in pp annihilation at rest into 7r + 7r~37r° |54| and 
pn annihilation at rest in deuterium into 27r + 37r~ []5S] (see also Given that mass 

(~ 1350 MeV) and width (~ 380 MeV) are compatible with / (1370), one might assume 
that / (1370) has been observed here in its pp decay mode. The large pp branching ratio 
points, however, to a large inelasticity in the nn S-wave around 1400 MeV. Thus / (1370) 
may split into two states, a QQ decaying to tttt, KK and rjr] and another state decaying 
to pp, possibly a molecule [ |57|j . 

The GAMS meson /o(1590) ]19[ decays to rjr)' and r\r\ with a relative branching ratio 
of 2.7 ± 0.8 |7J which, given the large error, is consistent with eqn. [T^ and [TJ5]for /o(1500) 
(see footnote 1). We shall assume that /o(1590) is either identical to /o(1500) or that it is 
the predicted ss state. We are therefore left with one nonet, two to three molecules and 
the supernumerary / (1500). 

If the ss member lies at 1600 MeV then from the linear mass formula 

^ = tf (61) 

3/o + ao - 4A 

where the symbols denote the particle masses, one obtains for the QQ nonet the singlet- 
octet mixing angle (59I7 2 ) (or 121°)). This is not far from ideal mixing, 9 = 35.3° (or 
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125.3°) . The relative decay rates of the ++ QQ mesons to two pseudoscalars are given in 
appendix A and hence a further consistency check on this proposed nonet may be applied 
when the partial decay widths eventually become known. Note that the experimental 
widths for the known QQ members of the proposed nonet (table [l|) are in good agreement 
with predictions, see appendix B. 

In the quark model, deviation from ideal mixing is due to mixing transitions between 
uu(dd) and ss which leads to a non-diagonal mass matrix in the flavour basis. The mixing 
energy (non-diagonal elements) is 

3A = f + f Q - 2K* (62) 

which vanishes for ideal mixing. For our ++ nonet one finds 37 MeV. The quark model 
also predicts the Schwinger sum rule 

(/o + /o) x (4^o - ao) " 3/o.fo " \( K * ~ «o) 2 = !j(a - K*f x cos 2 5 (63) 



where cos 2 <5 is the fractional QQ in the /q(1370) or /q("1600") wave function f58|. The 
fraction of glue in / (1370) is according to eqn. [JT} 

sm 2 5 = Y+~2{? < °' 33 - (64) 

Figure |TT] shows how the mixing angle 5 in eqn. ^ varies as the function of the ao mass 
for various /o(1370) masses. A consistent result is indeed obtained when m(/ (1370)) ~ 
1400 MeV and m(/ (1600)) = 1600 MeV. The fraction of glue in / (1600) depends on 
tuR 2 (eqn. 0). 



9 Conclusions 

We have argued that the properties of the / (1500) — / (1370) system are incompatible 
with them belonging to a quarkonium nonet. We suggest that /o(1500) is prominently a 
glueball mixed with the nn and ss quarkonia, the /o(1370) being dominantly Tin. 

We have shown that a reasonable scalar nonet can be built with the newly discovered 
ao(1450) setting the mass scale, and its width, together with that of the if* (1430), setting 
the scale of the nonet widths which are within 25 % of those expected by earlier quark 
model calculations. The width of /o(1500) is dramatically suppressed relative to these 
and the width of the /o(1370) may be also somewhat suppressed. These results are in 
line with these two states being partners in a glueball - QQ mixing scheme. It is this 
fortunate property that has enabled the /o(1500) to show up so prominently in several 
experiments where glueball channels are favoured. 

Further supporting these arguments we have shown that there appears to be no dra- 
matic intrinsic violation of flavour symmetry in decays involving gluons in a kinematic 
region where ++ or 2 ++ glueball bound states are expected to be negligible. The flavour 
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dependence of /o(1500) decays suggested that a significant mixing between G and QQ 
states is distorting the branching ratios in the ++ sector. We argued that the observed 
decay branching ratios could be due to the scalar glueball expected in this mass range, 
mixing with the two nearby QQ isoscalars, one lying below, the other above /o(1500). The 
partial decay widths of the lower state, /o(1370), are consistent with a mainly (uu + dd) 
state. Our hypothesis also implies that the (mainly) ss state lies in the 1600 - 1700 MeV 
region. 

The quantitative predictions of our analysis depend on the suppression of / (1500) 
decay to KK. Thus detailed study of pp — > ttKK can be seminal (i) in confirming the 
KK suppression, (ii) in confirming the K (1430) — > Kit and ao(1450) — > rjir and KK, (iii) 
in quantifying the signal for f (1370) and fo (1500) and (iv) in isolating the predicted ss 
member of the nonet. 

Clarifying the relationship between the Crystal Barrel /o(1500) and the /o(1590) of 
GAMS is important. Particular emphasis should be placed on the strength of the tttt 
branching ratio and the ratio of branching ratios r\r\j r n'n in light of its potentially direct 
significance as a test for glueballs. If the / (1550 ± 50) becomes accepted as a scalar 
glueball, consistent with the predictions of the lattice, then searches for the h and 
especially the 2 ++ at mass 2.22 ± 0.13 GeV |J may become seminal for establishing the 
lattice SIS db successful calculational laboratory. 
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Appendix A: 

Amplitudes for arbitrary pseudoscalar mixing angles 

The generalisations of the amplitudes for an arbitrary pseudoscalar mixing angle are as 
follows: 

Quarkonium decay 



(QQ\V\7T7r) 
{QQ\V\KK) 

(QQ\v\vv) 

(QQ\V\rrf) 




(65) 



for / = 



(QQ\V\KK) 
(QQ\V\irr)) 
(QQ\V\rrf) 




(66) 



for 1 = 1 and 



(QQ\V\Kn) 
(QQ\V\Kr)) 
(QQ\V\Krf) 



V3/2 

(psin0 — cos0/ y/2) / V2 
(pcos0 + sin0/ V2) I a/2 



(67) 



for I = 1/2. 



G QQQQ (graph §a) 



{G\V\nit) 
(G\V\KK) 

(G\V\rm) 
(G\V\rrf) 



R 

cos 2 + _R 2 sin 2 
cos0sin0(l — R 2 ) 



1 



(68) 
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G^GG (graph @b) 



(G|y|7T7r) = 
(G\V\KK) = 
(G\V\Tm) = (cos0V2 - Asin0) 2 
(Cj | V"| 7777') = (cos0\/2 — Asin0) (sin0\/2 + Acos0) 

(69) 



where 



in the pseudoscalar channel. 



_ (G(0-)\V\dd) 
X = (G(0-)\V\ss) m 



G^QQ (graph gc) 



(G|V|7T7T) = 1 

(G|y|KZ) = (p + uR 2 )/2 
(G\V\r]r}) = (cos 2 + upRhm 2 <f>) 
{G\V\rjrf) = cos0sin0(l — topR 2 ) 

(71) 
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Appendix B: 
Form Factors 



In the main body of the text we used rather simple forms for form factors, eqs (14). In 
order to test sensitivity to these assumptions we consider here the consequences of more 
structured form factors as arise when the dynamical effects of flux-tube breaking are 
included. The result is that momentum dependent multiplicative factors enter additional 
to those already in eqn (14); the general structure is discussed in table 2 and appendix 



B of ref |[44|| . In the approximation where the light hadron wavefunctions have the same 
scale parameter (the quantity (3 in eqn (14)) the structure of S— wave decay amplitudes 
(as for 3 P S + 1 S ) is 

q n 2g 2 g 2 , 



9p 2/ ^ 12/3 2 ' 

2 

The factor i n parenthesis, which was not present before, arises from the coupling 
of the 3 Pq of the initial QQ meson with the 3 Po of the qq pair created by the breaking 
of the flux-tube and which seed the decay. From detailed fits to meson spectroscopy and 
decays it is known that f3~ 2 ~ 5 to 6 GeV~ 2 and so the multiplicative factor (1 — J|s) 
does not significantly affect our analysis of the / (1500). 

In this appendix we have followed ref jPJ in imposing /3( 3 Pq) = 0.5 GeV and have used 



a more modern value J59i of (3 — 0.4 GeV for the 1 S'o and 3 P 2 states. Barnes et al |53j have 



made a similar analysis with (3 ~ 0.4 GeV as a preferred overall value and find similar 
results. The strategy is to use the known 3 P 2 decays to set the overall scale following the 



prescriptions in table 2 and appendix B of ref. f44| . 



(1) #(1370) 

Using r(/ 2 -> tttt) = 157 MeV yields, in MeV 

r(/ n -> tttt) = 270 ± 25 
KK = 195 ± 20 
r]T] = 95 ± 10 

where the r] is assumed to be a 50:50 mixture of ss and nn. The (1 — |C) factor has 
suppressed the 7T7t more markedly than the KK and rjrj, hence leading to a larger 7777 / 7T7T 
and KK/tttt ratio than in the main text. Nonetheless one sees that the r(/g l ) >> T(f£) 
still arises in line with the conclusion that a "narrow" /o width is out of line with a 
3 Po quarkonium state. This conclusion is reinforced by the expectation based on spin 
counting arguments that 7 2 (pp) > 7 2 (7T7t) and hence that there should be a non-negligible 
IX/q 1 — > p(7T7r) p ) in addition to the two-body channels. 
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(2) if *(1430) 

Using T(i^2 - ► Kir) ~ 50 MeV or r(/ 2 ) as above yields consistent similar results, namely 
in MeV 

T(K* -> Kir) = 200 ± 20 
Kr] ~ 
Kr]' = 15 ±20 

The Ifr/ has a large coupling for physical 7777' mixing angles but the width is very 
sensitive to phase space. We note that the Particle Data Group0] allow (7±10)% "non- 
i^7r" for the Kq decay and we have assigned this to Kr] 1 . 

These results should be compared with the experimental value T(K^) = 287±23 MeV. 



(3) a (1450) 

Using T(a2 — > KK) = 5.2±0.9 MeV (7] we obtain for the corresponding channel T(ao — > 
KK) ~ 110 MeV. Then with physical 7777' mixing angles we have r(?77r) ~ 90 MeV, 
T(r]'ir) ~ 80 MeV. There is considerable uncertainty in these widths, however, since if we 
were to normalise by T(/ 2 — > 7T7t) instead of by T(a 2 — > KK) we would find T(ao — > KK : 
7]7T : 7]'tt) ~ 200:160:140 MeV. Reflecting these uncertainties we can merely summarise by 

r(a -> 0-0~) theory = 390 ± 110 MeV 
r(a )e xp = 270 ±40 MeV. 

The general conclusion that T(J'q) > T(ao) > T(Kq) holds true. Furthermore we note 
that these results require that if the ao(1450) seen by Crystal Barrel is indeed 3 Pq (QQ), 
then comparable partial widths are expected for all of KK, rjir and t/tt. Quantifying these 
experimentally will be an important piece of the total strategy in clarifying the nature of 
these scalar mesons. 



(4) /<f(1.6-1.8) 

The (1 — ) can have a dramatic effect in the upper part of this range of masses. Using 
the r(/ 2 n — ► 7T7r) or r(/| KK) as normalisation, the mass dependence of the partial 
widths of a pure 3 Pq(ss) follows, in MeV 





1600 


1700 


1800 


r(/ * - kk) 


270 


155 


85 


rjrj 


45 


25 


20 


f]Vj 


195 


170 


190 



The effect of the momentum node is clearly seen in the KK and 7777 channels, whereas 
the 7777' maintains its strength. There is nothing in this pattern of widths that supports a 
3 Pq(QQ) interpretation of /o(1500). A detailed study of quarkonium decays with similar 



conclusions is in refp3]. 
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Figure Captions 



Figure 1: 7?- as a function of a for p = 1 (a) and p = 0.75 or 1.25 (b) for quarkonium 
decay (up to a common multiplicative factor). Dotted line: nn; dash-dotted line: KK; 
dashed line: 777/; solid line: 7777. 



Figure 2: x 2 f° r 2 ++ quarkonium decay as a function of p for various values of (5. The 5% 
C.L. limit is shown by the horizontal line. 



Figure 3: 7777 / 7T7T vs rjrj/KK invariant couplings per charge combination for various values 
of p < 1. The grey region where both ratios are larger than one is not accessible to QQ 
mesons unless ss production is enhanced. 
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Figure 4: ptan(ai) as a function of rjrj/irir (full curve) and 7777'/ 7T7T (dotted curve). The 
vertical arrows show the experimental ranges of Ri and i? 2 for / (1500) decay, left with 
form factor, right without form factor. 

Figure 5: p as a function of a. The full curves show the 5 and 10 % C.L. upper (+) and 
lower(-) limits dependence for the experimental value Ri, the dashed curves for R2. The 
dotted curves give the boundaries for the experimental value R3. The grey region shows 
the range allowed by the experimental data on / (1500) decay to -kit, rjr) and 7777'. The 
black region includes in addition the bubble chamber upper limit for KK. 

Figure 6: Contributions to gluonium decay: QQQQ (a), GG (b), QQ (c) and interpreta- 
tion as QQ mixing (d) involving the energy denominator E G — Eqq 

Figure 7: Predicted decay rates 7?- (up to a common multiplicative constant) (a) as a 
function of R 2 for QQQQ decay, (b) as a function of A for GG decay and (c) as a function 
of ujR 2 for QQ decays with p — 1. Solid line: irn; dashed line: KK; dash-dotted line: rjrj; 
dotted line: rjrj'. 

Figure 8: Glueballs, quarkonia and perturbations: (a) primitive QQ and (b) primitive 
glueball Go in flux tube simulation of lattice QCD; perturbation V\ (c) and V2 (d); the 
effect of V\ on QQ is shown in (e), and on G is shown in (f); the effect of V 2 on G is shown 
in (g). 

Figure 9: i? 3 versus 00 R 2 . 

Figure 10: Predicted ratio r as a function of 00R 2 for / (1500) decays. The solid lines show 
the four possible solutions from the Crystal Barrel results r\ and r 2 with the corresponding 
signs. The range allowed by the experimental errors is shown by the dashed lines for the 
"+ -" and "- -" solution. The solutions "+ -" and "- -" are compatible with uoR 2 in the 
range allowed by the KK suppression and radiative J/ip decay (area between the parallel 
lines, indicated by the grey flashes). A small breaking of chiral symmetry favors the "+ 
-" solution. 

Figure 11: Fractional contribution of QQ in the / (1370) wave function from the Schwinger 
sum rule as a function of a mass for various /o(1370) masses, assuming the (mainly) ss 
state to lie at 1600 MeV. The horizontal arrow shows the experimental uncertainty in the 
a mass and the vertical arrow shows the range allowed for |£| < 0.5. 
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